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O ■ Abstract 
O 

■ The nonmesonic weak decay of A hypernuclei is studied in a shell model 



framework. A complete strangeness-changing weak AA^ — > A^A^ transition po- 
tential, based on one boson exchange, is constructed by including the exchange 
of the pseudoscalar mesons vr, K, ry as well as the vector mesons p, iv, and K*, 



^ . whose weak coupling constants are obtained from soft meson theorems and 

SU(6)^. General expressions for nucleons in arbitrary shells are obtained. 
The transition matrix elements include realistic AN short-range correlations 
and NN final state interactions based on the Nijmegen baryon-baryon poten- 
tial. The decay rates are found to be especially sensitive to the inclusion of the 
strange mesons, K and K*, even though the role of kaon exchange is found to 
be reduced with recent couplings obtained from next-to-leading order Chiral 
Perturbation Theory. With the weak couplings used in this study the rates 
remain dominated by the pion-exchange mechanism since the contributions 
of heavier mesons either cancel each other or are suppressed by form factors 
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and short-range correlations. The total decay rate therefore remains in agree- 
ment with present measurements. However, the partial rates which are even 
more sensitive to the inclusion of heavier mesons cannot be reconciled with 
the data. The proton asymmetry changes by 50% once heavier mesons are 
included and agrees with the available data. 

I. INTRODUCTION 

In single A hypernuclei, a A hyperon can occupy any orbital in the hypernucleus since it 
is free from the Pauli exclusion principle due to its additional quantum number Strangeness. 
Hypernuclei are typically produced in some excited state through hadronic reactions such 
{K^,TT~) or {tt^,K'^) but can reach their ground state through electromagnetic gamma 
and/or nucleon emission. Eventually, they will decay through weak interaction processes 
which involve the emission of pions or nucleons but are nonleptonic in nature. 

A free A-hyperon has a lifetime of about 260 picoseconds and decays almost totally into a 
pion and a nucleon (A pn^ (~ 64%), A rnr^ (~ 36%)), with a release of kinetic energy 
of about 5 MeV to the nucleon along with a corresponding final momentum of about 100 
MeV/c. When the A is embedded in the nuclear medium, the phase space for the mesonic 
decay is greatly reduced since the A is bound by some 10 MeV for p-shell hypernuclei and 
up to 30 MeV in heavy hypernuclei. The final-state nucleon with its very low momentum 
thus becomes Pauli blocked, leading to a suppression of the mesonic rate by several orders 
of magnitude for heavy hypernuclei such as j^^Ph. Experimentally, however, one finds the 
lifetimes of hypernuclei to be roughly independent of A (see, e.g.. Fig. 14-1 in Ref. |jl|), 
though the data base is very poor, especially for systems with A >12. Therefore, the nuclear 
medium surrounding the bound A affects its weak decay by introducing new, nonmesonic 
decay modes, such as AA^ NN. Thus, hypernuclei larger than ^He decay mainly through 
these nonmesonic channels, where the A mass excess of 176 MeV is converted into kinetic 
energy of a final state of nucleons emerging with a momentum of about 400 MeV/c. As 
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a result, this process has a much larger phase space relative to the mesonic one, and the 
outgoing nucleons are not Pauli blocked. 

This was recognized more than forty years ago when it was suggested that the novel 
two-baryon decay mode AA^ —>■ NN could be understood in terms of the free-space decay 
mechanism A — * vrA^ with the exception that the pion now has to be considered virtual and 
is absorbed on a second nucleon bound in the hypernucleus. However, the large momentum 
transfer involved in the reaction leads to a mechanism that is sensitive to the short distance 
behavior of the amplitude and thus raises the possibility that the exchange of heavier mesons 
may play an important role. The production of these mesons would be below threshold 
for the free-space A decay, but they can contribute through virtual exchange in a two- 
baryon decay channel. As discussed further below it is the kinematic freedom of these 
additional boson exchanges that provides part of the motivation for this study. The fairly 
large momentum transfer also raises the hope that this reaction turns out to be insensitive 
to nuclear structure details and thus creates a suitable channel to investigate the weak decay 
mechanism. 

We would like to point out that there is another possible nonmesonic decay channel, the 
two-nucleon induced process ANN NNN, where the virtual pion emitted at the weak 
vertex is absorbed by a pair of nucleons which are correlated through the strong force. This 
mechanism was first investigated in Ref. where it was suggested that its magnitude could 
be comparable to AA^ NN. However, a reanalysis with more realistic assumptions [^,^ 
reduced its contribution to 10-15% of the total nonmesonic decay rate. Its relevance lies 
mainly in its potential to renormalize the experimentally measured partial ratios; we will 
return to this point later when discussing results. 

The main goal in studying the weak nonmesonic decay channel is to gain insight into 
the fundamental aspects of the four-fermion, strangeness changing weak interaction. Most 
of the earlier work on nonleptonic weak processes proceeded directly from some model weak 
Hamiltonian and computed experimental observables that could then be compared to mea- 
surements. More recently the approach has been divided into a two-stage process. The first 
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step starts with the Standard Model electroweak Hamiltonian at the sub-hadronic level and 
takes into account QCD corrections at short distances, yielding a so-called effective weak 
Hamiltonian. In relation to the mass of the W-boson of 80 GeV, these would be low momen- 
tum transfer processes leading to a zero-range or contact interaction. The Cabbibo theory 
combined with strong-interaction corrections would result in V — A weak interaction and 
presumably predict the relative strength of the AS'=0 and AS'=1 transition. Thus, hadronic 
weak matrix elements of the form {MB'\Hw \ B) could eventually be calculated. The second 
stage involves using these weak vertices as a starting point for effective nuclear two-body 
operators that are then implanted into the nucleus with the usual nuclear many-body wave 
functions. It is the latter task which is the main subject of this paper. 

The most important information regarding the nonleptonic weak interactions of hyperons 
comes from their free-space decays. However, while these processes have been well-measured 
and understood phenomenologically, a more basic understanding in terms of the underlying 
degrees of freedom is still lacking. The situation is similar for the nonmesonic decay modes, 
except that the experimental data are more unsatisfactory. The early measurements were 
based on bubble chamber experiments and emulsion works, who suffered from low precision, 
poor statistics and difficulties with the identification of the particular hypernucleus, except 
for the very light hyperfragments. These experiments were nevertheless able to establish 
the first limits on hypernuclear lifetimes, albeit with large uncertainties. Their techniques 
were mostly aimed at measuring the emission of the vr" from the mesonic decay channel; 
the analyses of the nonmesonic decay modes were greatly inhibited by the low hypernuclear 
production rates, the limited spacial resolution of early detectors and the presence of one 
or more neutral particles (especially neutrons) in the final state. In recent years a series of 
counter experiments carried out at BNL and KEK improved the quality of data on the non- 
mesonic decay modes using pion and kaon beams. In contrast to these measurements LEAR 
at CERN explored hypernuclei in the A=200 mass region following antiproton annihilation 
1^. The heavy hypernuclei were identified through their delayed fission events which were 
attributed to the nonmesonic decays. While this interpretation is fraud with difficulties it 
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resulted in lifetimes similar to the those of p-shell hypernuclei, thus raising the intruiging 
possibility that the AA^ — > A^A^ mechanism saturates already in the A=12 system. On the 
other hand, employing the {K^,tt^) production mechanism with direct timing techniqes, 
the Brookhaven AGS measurements resulted in better data on the hypernuclear lifetimes 
and branching ratios of mesonic as well as nonmesonic decay channels of ^^C, ^^B and the 
different helium hypernuclei The group at KEK used the (vr^, K^) production reaction 
to not only measure total and partial decay rates 0, but also to use the induced polarization 
in j^C to determine for the first time the difference in the number of protons emitted along 



the axis of polarization compared to the number ejected in the opposite direction |jTO[. This 
asymmetry is a direct consequence of the presence of parity-violating components in the 
weak decay mechanism. 

Parity violation in hadronic systems represents a unique tool to study aspects of the 
nonleptonic weak interaction between hadrons. The nonmesonic process resembles the weak 
AS=0 nucleon-nucleon interaction that has been explored experimentally in parity-violating 
NN scattering measurements by measuring the asymmetry of longitudinally polarized pro- 
tons. However, the AA^ — > A^A^ two-body decay mode contains more information since it can 
explore both the parity-conserving (PC) and the parity-violating (PV) sector of the AS=1 
weak baryon-baryon interaction while in the weak NN system the strong force masks the 
signal of the weak PC interaction. On the other hand, the free process AA^ —>■ NN cannot 
be accessed experimentally which complicates the interpretation of the nonmesonic decay 
rates since the reaction mechanism has to be studied in the environment of hypernuclear 
structure. 

A number of theoretical approaches to the AA^ NN decay mode have been developed 



over the last thirty years which are more extensively reviewed in Ref. [|Tl|. The early phe- 
nomeno logical analyses by Dalitz et. al. |jl2[| provided the general nonrelativistic structure of 
the AA^ — > A^A^ amplitude which was then related to decay rates of s-shell hypernuclei using 
certain simplifying assumptions. The AS* = weak nucleon-nucleon interaction at low and 
intermediate energies has generally been described in a meson exchange model involving one 



strong interaction vertex and one weak one; the same basic assumption has been used for 
a microscopic description of the AS=1 AN NN mechanism. Early calculations based 



on the one-pion exchange (OPE) were due to Adams [0, modifications of the OPE due to 



strong interactions in the nuclear medium were suggested in Ref. to account for many- 
body nuclear structure effects. At the very least, the OPE mechanism can be expected to 
adequately descibe the long-range part of the AA^ — A^A^ interaction. The first attempts 
to include heavier bosons, at first the p meson, — again in complete analogy to the AS'=0 
NN interaction — were presented in Ref. ||l5| , p!6| . There were several conference papers by 



Dubach et al ^Tj showing results of preliminary calculations with a full meson exchange 



potential; a more detailed account of their calculations has recently become available [jTH 
Finally, for completeness we mention that a study of AA^ NN based on quark rather than 
meson degrees was carried out in Ref. Their model separates the process into a long- 



range region, to be described by OPE, and a short-range region, modelled by a six-quark 
interaction with suitably adjusted parameters. This idea was revived in Ref. ||21[]. Results 
for light hypernuclei showed that the total rates preferred a relative plus sign between quark 
and OPE contributions, while the neutron to proton ratio was better reproduced with a mi- 
nus sign. In the calculation of Ref. both OPE and OKE (one-kaon exchange) amplitudes 
were considered, in addition the the quark ones, and significant cancellations between OPE 
and OKE contributions were found. Finally, there have been some attempts to incorporate 
the exchange of the a and p mesons from the point of view of a correlated two-pion exchange 
23| , |2^ , where the weak vertex was obtained through the coupling of the two pions to the 



a oi p and intermediate N and S baryon states. Both calculations find an enhancement of 
the central transitions leading to a moderate increase of the neutron to proton ratio. 

The motivation for this work is two-fold. First, in contrast to most previous investigations 
performed in nuclear matter this study analyses the nonmesonic hypernuclear decay in a shell 
model framework. Spectroscopic factors are employed to describe the initial hypernuclear 
and final nuclear structure as well as possible. Rather than constraining the AN system 
to only /=0, all possible initial and final relative orbital angular momenta are included. 
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To reduce the uncertainties regarding initial and final short-range correlations (SRC) we 
use realistic AN and NN interactions based on the Nijmegen baryon-baryon potential. The 
nuclear structure details are thus treated with as few approximations and ambiguities as 
possible. We emphasize here that our current treatment is nonrelativistic, in contrast to 
our previous works pB| , |2B| . As discussed in detail in Refs. P7| , p8| , one has to pay special 
attention to the effects of short-range correlations in a relativistic treatment of two-body 
matrix elements. Secondly, our calculations are performed in a full one-boson-exchange 
model that includes not only the long-ranged pion but also contributions from the other 
pseudoscalar mesons, the rj and K, as well as the vector mesons p,uj and K*. Since nuclear 
structure uncertainties have been eliminated or are minimal, we can use our framework 
to draw conclusions regarding the sensitivity to the underlying weak baryon-baryon-meson 
couplings. The many-body matrix elements and two-body amplitudes are evaluated in 
Section II and III respectively. The form of the meson exchange potential is presented in 
Section IV for the pseudoscalar and vector mesons. Section V discusses our calculation of 
the coupling constants. Our results are given in Section VI where we discuss the influence 
of each meson on the total and partial rates and the proton asymmetry. Our conclusions 
are presented in Section VII. 



II. DECAY RATE AND ASYMMETRY 

The nonmesonic AA^ —>■ NN decay rate is given by 

r_ = / ^ / (^E(2-)^(M^ -En-E,- E,) \ M p (1) 

where = {i/jr; Pk S Ms T Mt\ 0\n^nn \aA) is the amplitude for the transition from an 
initial hj^ernuclear state to a final state which is divided into a two-nucleon state and a 
residual {A — 2)-particle state. The quantities Mh, Er, Ei and E2 are the mass of the 
hypernucleus, the energy of the residual {A — 2)-particle system, and the total asymptotic 
energies of the emitted nucleons, respectively. A transformation to the center of mass (P) 
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and relative momentum (k) of the two outgoing nucleons is already implied in Eq. (|l|). The 
sum J2 indicates an average over the initial hypernucleus spin projections, Mj, and sum over 
all quantum numbers of the residual {A — 2)-particle system, as well as the spin and isospin 
projections of the exiting nucleons. We follow Ref. by assuming a weak coupling scheme 
where the isoscalar A in an orbit = {tt-a, l\, s\,j\, m\} couples to only the ground state 
wave function of the nuclear {A — 1) core 

= E {Ja^aJcMc I JiMi) I K/ASA)jAmA) I Jc'McTjTs^) . (2) 

TTlA Mc 

Employing the technique of coefficients of fractional parentage, the core wave function is fur- 
ther decomposed into a set of states where the nucleon in an orbit ajy = {un, 1^, sjy, j^, tjin} 
is coupled to a residual {A — 2)-particle state 

I JcMcTiT^j) = E {JcTiil JRTR,jNtN)[ I Jr, Tr)x \ {nNlNSN)jN, tN)YTiTl'^ 

= E i-JcTiil JRTR,jNtN) 
JrTrJn 

X E E {JRMRjNfn^ I JcMc){TnT3^tr,h^ \ Ti^^) 

Mr mjv Ts^ is, 

X I JrMr) I TrTsj^) I {nNlNSN)jNmN) I tNts,) , (3) 

where t]\r = 1/2. The spectroscopic factors, S = {A — l){JcTf{\ JRTji,i]^tiq)'^, appropriate 
for the decay of a^C, are taken from Ref. [^] and are listed in Table 4 of Ref. |26|. Taking 
into account that the A decays from a Za = state and working in a coupled two-body spin 
and isospin basis, the nonmesonic decay rate in Eq. can be written as 

^nm ~l~ ; (4) 



where 



X E I {TrTsrUs^ I T.Ts,) r 
Mi ^ 



X 



TTz ^ ^ mt^Mc jN 



1 11 

X E {JRMnjNmN \ JcMc) E {Jn^n \ iNmn^-rris,^) E ('S'oMso | -^a -m^^^) 

2 

X ^(ToTaji-i ^tg.) tA7v-.7V7v(^, M5, T, M^, ^o, M^^, To, Tg^, P, k) , (5) 
T0T30 ^ ^ 

with ^3- = 1/2, ti = —1/2, t2 = 1/2 for the p-induced rate {Ap np) and ^3. = —1/2, 
ti = —1/2, ^2 = ~l/2 for the n-induced rate (An nn). Equation @ is written in terms 
of the elementary amphtude tAAf^TVAf, which accounts for the transition from an initial AN 
state with spin (isospin) Sq (Tq) to a final antisymmetric NN state with spin (isospin) S (T). 
The details on how this two-body amplitude is calculated are given in the next section. Note 
that the A has been assumed to act as a | 1/2 — 1/2) isospin state which is coupled to the 
nucleon to total isospin Tq. As explained in the next section, this is the way to incorporate 
the change in isospin A/ = 1/2 induced by the weak transition operator. 

At the kinematic conditions of the {'n'~^,K~^) reaction carried out at KEK, the hypernu- 
cleus is created with a substantial amount of polarization in the ground state. Due to the 
interference between the parity conserving and parity violating amplitudes, the distribution 
of the emitted protons in the weak decay displays an angular asymmetry with respect to 
the polarization axis given by 

cx{x)=ao{l + PyA,{x)) (6) 

where Py is the hypernuclear polarization, created in the production reaction, such as 
(tt^, K^) at KEK and BNL or (7, K^) at CEBAF ^0|, and the expression for the asymmetry 

AM- 3 TriMSyM^) _ 3 . 

^^^^^"j+1 Tr{MM^) - J+1 E^(M,) 

Mi 



shown, for instance, in Ref. defines the asymmetry parameter, Ap, characteristic of 
the hypernuclear weak decay process. The asymmetry in the distribution of protons is thus 
determined by the product PyAp. In the weak coupling scheme, simple angular momentum 
algebra relations relate the hypernuclear polarization to the A polarization 



-Py if J — JC — I 



PA={ ' ' ^ , (8) 

Py if J=Jc + l 

where Jc is the spin of the nuclear core. It is convenient to introduce the intrinsic A 
asymmetry parameter 

f-^A, if J = Jc-i 

OA = < ^ , (9) 

[ Ap if J = Jc + 1 

such that PyAp = p^a^, which is then characteristic of the elementary A decay process, 
AA^ — s> A^A^, taking place in the nuclear medium. 



III. TWO-BODY AMPLITUDES 

In this section we describe how we evaluate the elementary two-body transition amplitude 
tATV^TVAT which, as shown in Eq. (|), contains the dynamics of the weak decay process. In 
the first place, it is necessary to rewrite the product of two single particle wave functions, 
(ri I a\) and (r2 | a^), in terms of relative and center-of-mass coordinates, r and R. In 
the present work, the single particle A and N orbits are taken to be solutions of harmonic 
oscillator mean field potentials with parameters 6a = 1.87 fm and = 1.64 fm, respectively, 
that have been adjusted to experimental separation energies and the ^^C charge form factor. 
Assuming an average size parameter b = {bx + bN)/2, using Moshinsky brackets and working 
in the LS representation, the product of the two harmonic oscillator single particle states, 
^nimi'^i) <l>^;/„/(r2), can be transformed to a linear combination of products of relative 
and center-of-mass wave functions. Since the A is in a /a = shell, we obtain 

(t) *™ (f ) - *™ «■ (^) • 

when the nucleon is in the s-shell and 
<o (^) ) ^ {<. (^) (^) - (^) ^T. (^) ) 
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when the nucleon is in the p-shell. 

As for the final NN state, the antisymmetric state of two independently moving nucleons 
with center-of-mass momentum P and relative momentum k reads 

{Rr\Pk S Ms TMt) = -^e^^^ (e^'^^ - (-l)^+^e-'^r) ■ (12) 

To incorporate the effects of the NN interaction a substitution of the plane wave by a 
distorted wave 

needs to be done. Therefore, the amplitude Ian^nn of Eq. (H) can be decomposed in terms 
of amplitudes which depend on CM. and relative quantum numbers 

X{NrLrNnLn,hlN)t^htNN''' (14) 

NrLrNRLR 

where X{NrLrNRLR,l\lN) are the Moshinsky brackets which for 1\=In=0 are just 
X(l 1 0, 0) = 1, and for = 1 are X{1 110,01) = 1/^2 and X{1 101,01) = -1/V2, 
as can be seen from the decomposition of the wave function above. The matrix elements 
tA^N^NN^'' are given by 

fNrLrNnLn _ J_ [ j3 [ ,3 ~~iPn^*(^^ iS ]Ty (-.^CM ( ^ \ ^rcl ( ^So To 

tAN^NN - "^y "-re y^kirjXMsXTa ^ [^) "^NuLr ^rLr y^i^j ^Mso^ng 

= (^V) ^rel , (15) 



V2. 



with 



where, for simplicity, we have only shown the direct amplitude corresponding to the first 
term of Eq. (|1^). The function ^%^Lii ^"^"^^ Fourier transform of the AN center-of- 

mass wave function and trei is the expectation value of the transition potential ^(r) between 
AN and NN relative wave functions. In the next section it is shown how the potential V^(r) 
can be decomposed as 
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nr)=EEKl'Hr)Oj«, (17) 



i a 



where the index i runs over the different mesons exchanged and a over the different spin 
operators, Oa G (1, cricr2, 5'i2(r) = 3crircr2r — cri(T2, cr2r, [cti x (T2]r), which occur in the 
potential V{r). The isospin operator, depends on the meson and can be either 1 for 
isoscalar mesons (?7, lu), tiT2 for isovector mesons (vr, p) or a hnear combination of 1 and 
T1T2, with the coefficients depending on the particular spin structure piece of the potential, 
for isodoublet mesons (K, K*). The radial parts V^'^ (r) are discussed in the next section. 

By performing a partial wave expansion of the final two-nucleon wave function and 
working in the (LS*) J-coupling scheme, the relative AA^ NN amplitude, Uei, can be 
further decomposed 

= ^EE E ^7^^-''\LMLSMs\JMJ)YLMM{LrML^SoMs,\JMJ) 

v2 i a LL'J 

X {{L'S)JMj I 6, I iLrSo)JMj){TT, I f ^ I TO^) J r'drni'iK,r)V^'\rWNlL^r) , (: 

where ^^^rLri''^) stands for the radial piece of the H.O. wave function. The explicit expressions 
for the expectation value of the spin-space piece, {{L'S)JMj \ Oa \ {LrSQ)JMj), can be 
found in the Appendix. 

The function f^]^,{kr,r) is the scattering solution of two nucleons moving under the 
influence of the strong interaction, for which we consider the updated version of the Reid- 
soft core potential |^T|, given in Ref. H], and the Nijmegen H] NN potential. Such a wave 



function is obtained by solving a T-matrix equation in momentum space and in partial wave 
decomposition following the method described in Ref. |3^. The tensor component of the 
NN interaction couples relative orbital states {L and L') having the same parity and total 
angular momentum as, for instance, the ^Si and ^Di channels. Therefore, starting from an 
initial Lr orbital momentum, the weak transition potential produces a transition to a L' 
value, which mixes, through the subsequent action of the strong interaction, with another 
value of orbital angular momentum, L. In Table | we present all the possible final states 
starting from initial AN states having either L,, = or 1 and for the central (AS = 0, 
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AL = 0), tensor (AS = 2, AL = 2) and parity violating (AS = 1, AL = 1) pieces in which 
the transition potential can be decomposed. 

In the absence of final state interactions (FSI), the NN wave function in Eq. (|T8|) would 
reduce to a spherical Bessel function 

^LL'{kr,r) = 6L,L'jL{krr) . (19) 

We note that the procedure followed here to include FSI between the two emitted nucleons 
differs from our previous works P^p5[, where the non-interacting NN pair, represented by 



a Bessel function in the final state, was multiplied by an average NN correlation function 

/FSi(r) = l-Jo(gcr), (20) 

with Qc = 3.93 fm~^, which provides a good description of nucleon pairs in ^He as 
calculated with the Reid-soft core interaction . 



To account for the AN correlations, which are absent in the independent particle model, 
one should replace the harmonic oscillator AN wave function, $5^'^^(r), by a correlated AN 
wave function that contains the effect of the strong AN interaction. Such wave functions were 
obtained from a microscopic finite-nucleus G-matrix calculation using the soft-core and 



hard-core Nijmegen models [38|. In Ref. [27| we showed that multiplying the uncorrelated 



harmonic oscillator AN wave function with the spin-independent correlation function 

/(r) = (l - e-'^'/'^y + fer^e-^'/'^' , (21) 

with a = 0.5, b = 0.25, c = 1.28, n = 2, yielded decay rates slightly larger than those 
obtained with the numerical Nijmegen soft-core correlations but slightly smaller than those 
computed with the Nijmegen hard-core potential. Since the deviations were at most 10% the 
above parametrization can be used as a good approximation to the full correlation function. 



IV. THE MESON EXCHANGE POTENTIAL 



The transition AA^ — NN is assumed to proceed via the exchange of virtual mesons 
belonging to the ground-state pseudoscalar and vector meson octets. As displayed in Figs. 
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|T|(a) and (b), the transition amplitude involves a strong and a weak vertex, the later being 
denoted by a hatched circle. 

A. Pseudoscalar mesons 

While there exist several strong meson-exchange potentials which, through fits to NN 
scattering data, provide information on the different strong NN-meson vertices, only the pion 
vertex is known experimentally in the weak sector. The weak hamiltonian is parametrized 
in the form 

= iGprnlMA^ + B^^,)Tcf>-^Pj, (?) , (22) 

where Gprn'^ = 2.21 x 10~^ is the weak coupling constant. The empirical constants = 1.05 
and 5^ = —7.15, adjusted to the observables of the free A decay, determine the strength 
of the parity violating and parity conserving amplitudes, respectively. The nucleon, lambda 
and pion fields are given by ■j/'n, "^Aa and (f)'^, respectively, while the isospin spurion (?) is 
included to enforce the empirical A/ = 1/2 rule observed in the decay of a free A. 
For the strong vertex, we take the usual pseudoscalar coupling 

nlj,^ = igNNTTipNl^rcp^'ipN , (23) 

which is equivalent to the pseudovector coupling when free spinors are used in the evalua- 
tion of the transition amplitude. The nonrelativistic reduction of the free space Feynman 
amplitude is then associated with the transition potential. In momentum space, one obtains 

V;N^-G.™^^(i + ^..<,)^ (24) 

where q is the momentum carried by the pion directed towards the strong vertex, g = ^f^N^ 
the strong coupling constant for the NNvr vertex, fi the pion mass, M the nucleon mass and 
M the average between the nucleon and A masses. The operators A and B, which contain 
the isospin dependence of the potential, read 
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A = TiTa 

B = B^ T^T2 (25) 

We note that the nonrelativistic approach of the present work differs from our previous 
works [p5| , p6[ , which were based on a relativistic formahsm. It was found that the supression 



of the matrix elements due to short range correlations was larger by about a factor of two 
to what was obtained in standard nonrelativistic calculations |]T^,|T3|,[T^,|TB[ . In Ref. [|27| it 



was shown that, if one uses the same nonrelativistic correlation function, the relativistic 
and nonrelativistic schemes were not giving the same correlated potential obtained through 
the standard nonrelativistic reduction. In the relativistic approach, the correlation function 
was applied to the Feynman amplitude before the nonrelativistic reduction was carried 
out, whereas in the nonrelativistic procedure the correlation function was applied after 
the reduction of the free Feynman amplitude was obtained. The difference between the 



two methods was studied in Ref. |]28|, where it was shown explicitly that the relativistic 
framework together with a standard nonrelativistic correlation function lead to additional 
contributions in the correlated transition potential which produced the larger supression of 
the decay rates reported in Refs. [25,2^. For lack of relativistic correlation functions we 
adopt a nonrelativistic formalism in this paper. 

The other mesons of the pseudoscalar octet are the isosinglet eta (r/) and the isodoublet 
kaon (K). The strong and weak vertices for these mesons are 

^NN,? = ^^nn,,^n750''^n (26) 
nZ, = tGFm'ji^ (A, + 5,75)0>A (?) (27) 

^ANK = 'i'9ANKi>Nl5 (l)^1pA (28) 

nZ^ = iGpml (?) (Cr + (0^) Vn + V^nV^n {Dl- + (0^)^ (?)] (29) 

where the weak coupling constants cannot be derived from experiment. In the present work 
we adopt the approach of Refs. [0,^], presented in the next section. 

We note that the isospurion (?) appearing in the former equations is used to enforce the 
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empirical AJ = 1/2 rule. The particular structure of the K weak couplings reproduce the 
vertices shown in Fig. ^ 

The corresponding nonrelativistic potentials for the transition AA^ —>■ NN are analogous 
to Eq. (^) but making the following replacements 

9 5'nN77 

fj. rriri 
A > Ari 



B^B^, (30) 



in the case of ry-exchange, and 



9 ~^ A'ank 

/i — * mK 
A - 



2 




2 


2 




/^PC 

- 2 



M 

B^[^ + Dl^ + ^T,T,] , (31) 

in the case of K-exchange, where the factor = corrects for the fact that the nonrelativistic 
reduction of the strong ANK vertex gives a factor 1/M instead of 1/M. Performing a Fourier 
transform of the general expression given in Eq. (^41) and introducing the tensor operator 
'S'i2(r) = 3crifcr2r — cricr2, it is easy to obtain the corresponding transition potential in 
coordinate space, which can be divided into central, tensor and parity violating pieces. The 
explicit expressions are given at the end of this section. 

B. Vector mesons 

A number of theoretical studies in recent years have investigated the contribution of the 



p-meson to the AA^ — > A^A^ process [|T^,|T^,^]. The weak ANp and strong NNp vertices are 
given by |jl5 
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= GFml (a,j^ - p/-^ + e,ri^) rp^i,^ (?) (32) 

^Lp = ?N {gUl' + ^i^^"'?-) ^P.V'N , (33) 

respectively, where the four momentum transfer q is directed towards the strong vertex. The 
values of the strong and weak coupling constants are given in the next section. 

The nonrelativistic reduction of the Feynman amplitude gives the following p-meson 
transition potential 

n (a + /3)(Fi + F2) , 



, .£(Fi + F2) , , \ 1 



with \i = nip, Fi = g^^p, F2 = S'Jnp the operators a, [3 and e 

a = Up T1T2 
(3 = (3p T1T2 

e = Ep T1T2 (35) 

contain the isospin structure. Using the relation (cri x q)(cr2 x q) = (cricr2) q^ — (criq)(cr2q) 
and performing a Fourier transform of V^(q), one obtains the corresponding transition po- 
tential in coordinate space, which, as in vr exchange, can be divided into central, tensor 
and parity-violating pieces. Furthermore, the p-meson central potential can be further de- 
composed into a spin independent and a spin dependent part [0. Due to the different 
models employed for the weak ANp- vertex [ll5| , p!6| , |39| , different calculations have yielded 
widely varying results. However, all studies until now have only included the tensor piece 
of the parity-conserving p-exchange term motivated, in part, by the observation that this 
is the most important contribution to the vr-exchange potential. We recently demonstrated 



3^ that the central piece of the p-exchange is in fact larger than its tensor interaction, an 
observation that can be traced to the fact that the p-exchange diagram has a much shorter 
range than the vr-exchange potential. It is therefore important to explicitly keep all pieces 
of the potential for the vector mesons. 
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The other vector mesons considered in this work are the isoscalar uj and the isodoublet 
K*, for which the weak and strong vertices can be written as 



n 



w 



n 



w 



+ 



Ck^^n (?) (0f )^ 7^75^N + I^^y?N7^5^N(0f )^ (?); 



(36) 
(37) 
(38) 



+ 



(39) 



Note that the K* weak vertex has the same structure as the K vertex, the only difference 
being the parity conserving contribution which has two terms, related to the vector and 
tensor couplings. The nonrelativistic potential can be obtained from the general expression 
given in Eq. (|33) making the following replacements 



/i 

a 



m„ 



V 



T 

9nnuj 



a. 



in the case of u exchange, and 



(40) 



Fi 
F2 
a 



V 

5'ank* 



T 

S'ank* 



^PC,V ^PC, 



-TIT2 
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/OPC,T ^PC,T 

P ^ ^K* H ^^1^2 

/(^pv ^pv \ f,/r 

' ^' ' Z^- + %rir, i, (41) 



2 ^ 2 M 



for the K* meson exchange. 



C. General Form of the potential 

The Fourier transform of the general Eqs. (24) and (34) leads to a potential in configu- 
ration space which can be cast into the form 

V{r) = EEK.«(r) = EEK!'Hr)o./« 

i a i a 

= E \V^P{r)W + V^{r)cT,aM + V}f> (r) 3^(^)1? + 

i 

+ (nV2 . r + (1 - [cr^ x a^] ■ r) F#(r)/lj|.] , (42) 

where the index i runs over the different mesons exchanged [i = 1,...,6 represents 
n, p,K,K* ,ri,uj) and a over the different spin operators denoted by C (central spin inde- 
pendent), SS (central spin dependent), T (tensor) and PV (parity violating). In the above 
expression, particle 1 refers to the A and n* = 1(0) for pseudoscalar (vector) mesons. In 
the case of isovector mesons (vr, p) the isospin factor is T1T2 and for isoscalar mesons {ri,uj) 
this factor is just 1 for all spin structure pieces of the potential. In the case of isodoublet 
mesons (K, K*) there are contributions proportional to 1 and to T1T2 that depend on the 
coupling constants and, therefore, on the spin structure piece of the potential denoted by a. 
For K-exchange we have 



f(3) 


= 






f(3) 
^SS 


f(3) 

— Irp 


2 




f(3) 


2 







T1T2 

(43) 



and for K* exchange 
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f(6) 

r(6) 
'SS 

f(6) 



+ D 



r(6) 



+ 



a': 



PC,V 



2 " 2 



-TIT2 



PC,T^, 



(44) 



The different pieces V^^'> , with a = C, S'S', T, , are given by 



6{r) 



K^lVss{r,fii) 



Q-I^i^ / 1 \ r\ 



- 1 + — ^K';>yVpvir,fi, 
Anr \ fiir J 



(45) 
(46) 

(47) 

(48) 



where /i, denotes the mass of the different mesons. It is these expressions that are inserted 
in Eq. ([T8|) to compute trei numerically. The expressions for Kj^^ , which contain factors and 
coupling constants, are given in Table |I[ 

A monopole form factor -Fi(q^) = (A^ — /x^)/(A^ + q^) is used at each vertex, where 
the value of the cut-off, Aj, depends on the meson. We take the values of the Jiilich YN 
interaction displayed in Table ^ of Section V, since the Nijmegen model distinguishes 
form factors only in terms of the transition channel. The use of form factors leads to the 
following regularization for each meson 

2 



Voir; /X,) - Vc{r; /x.) - Vc{r; A,) - A.^^^^^^ (l 
Vssir; /i.) ^ Vssir; /x,) - Vss{r- A,) - A^^^l^J^^ ~ 



Vrir; Vrir; - Vrir; Ai) - A, 



Vpv{r; /ii) Vpv{r; /ij) - Vpy(r; Ai) 



V^^e^^ / 1 
2 47r V Air 

2 Air 



(49) 
(50) 
(51) 
(52) 



where V^(r; Aj) has the same structure as V^(r; /Xj), defined in Eqs. (|45|)-(^8D, but replacing 
the meson mass jUj by the corresponding cutoff mass Aj. 
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V. THE WEAK COUPLING CONSTANTS 



The starting point for describing the weak decay of strange particles has been the fun- 
damental Cabbibo Hamiltonian based on the Current Current assumption 

G 



H = -^J d^xr{x)Jl{x) + h.c. (53) 



with 



+ M(x)7a(l - 75) cos^c + ■s(x) sin^c-) , (54) 

where Oq is the Cabbibo angle, Gp the weak coupling constant and we take the Bjorken 
and Drell convention for the definition of 75 As is well known, terms proportional 

to cos describe, for instance, the neutron /5-decay while the contributions proportional 
to sin6'c lead to the semileptonic decay of hyperons and kaons. The AS'=1 nonleptonic 
decays are governed by terms proportional to sin dc cos dc which consist of products of a 
current between u and d quarks (A/=l) and a current betwen u and s quarks (A/ = 1/2). 
Thus, since terms in sin ^c* cos 6'c describe transitions with A/ = 1/2 and 3/2 with equal 
probability, the empirical A/ = 1/2 rule indicates the presence of some dynamical effect 
related to QCD corrections that suppresses the AJ=3/2 components of the Hamiltonian. 

In order to obtain hadronic weak matrix elements of the kind {MB'\Hu)\B) , where M 
can stand for pseudoscalar or vector mesons and B for baryons, it has been convenient to 
express the effective weak Hamiltonian in terms of the SU(6)^ symmetry that unites them. 

The AS=1 weak nonleptonic Hamiltonian can be written in SU(3) tensor notation: 



Hu, = ^ cos 9c sin 9c{Jl„ J^} + h-c (55) 



where J^^ = (V^ — A^Yj is the weak hadronic current with SU(3) indices i and j. As shown 
in Ref. [|1^,^,^ the weak vector and axial currents can be expressed in terms of SU(6)^ 
currents. Since the Hamiltonian is the product of two currents, each belonging to the 35 
representation, one can expand 
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35 35 = 1, ® 35s © 189, © 405, © 35„ © 280„ © 280^ , (56) 

which allows extraction of the parity violating (PV) and parity conserving (PC) pieces of 
the Hamiltonian: 

Hpc ■■ h ® 35, © 189, © 405, (57) 

Hpy : 35a ® 280„ © 280^ . (58) 

Each of the possible ways of coupling baryons to mesons within the SU(6)^ symmetry intro- 
duces a reduced matrix element that can either be fitted to experimental data or calculated 
microscopically from quark models. Below we discuss the PV and PC amplitudes separately. 

A. The PV amplitudes 

The traditional approximation employed to obtain the PV amplitudes for the nonleptonic 
decays B ^ B' + M has been the use of the soft-meson reduction theorem: 

WMB' M,{q)\H^,\B) = -±{B'\[F!, H,,]\B) = -±{B'm, H^,]\B) (59) 

where q is the momentum of the meson and Fi is an SU(3) generator whose action on a 
baryon Bj gives: 

Fi\Bj) = ifijk\Bk) . (60) 

Since the weak Hamiltonian is assumed to transform like the sixth component of an 
octet, a term hke {Bk\H^\Bj) can be expressed as: 

{Bk\Hl\Bj) = iFfejk + Ddejk (61) 

where fijk and dijk are the SU(3) coefficients and F and D the reduced matrix elements. 

With the use of these soft-meson techniques and the SU(3) symmetry one can now relate 
the physical amplitudes of the nonleptonic hyperon decays into a pion plus a nucleon or a 
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hyperon, B ^ B' + ir, with the unphysical amphtudes of the other members of the meson 



octet, the kaon and the eta. One obtains relations such as [H2|: 



2- ^2 



{nv\H^.\A) = ^h'_ (62) 

where Eq" stands for (p7r°|/7pv|S^), the PV amphtude of the decay E"*" pn^, which is 
experimentally accessible. We have used the standard notation according to which the 
hyperon and meson charges appear as superscript and subscript, respectively. We point 
out that using the isospin structure of the potential defined in the previous chapter, the 
{NK\Hp^\N) matrix elements are connected to the coupling constants C^^ and of Eq. 
(0) via 

{nK^\H^,\p) = Cl^ (63) 

As shown above, the symmetry of SU(3) allows connecting the amplitudes of the physical 
pionic decays with those of the unphysical decays involving etas and kaons. SU(6)^, on 
the other hand, furthermore permits relating the amplitudes involving pseudoscalar mesons 
with those of the vector mesons. For details we refer the reader to Refs. we just list 

here the final relations in terms of the coupling constants defined in the previous section, 
rather than matrix elements 
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^K* = + ^a^ (64) 

The numerical values of the constants are given in Table |T|. Note, that an additional 
parameter, ot, is present in the coupling constants for the vector mesons. This coupling, 
which is very small in the case of pion emission due to PCAC, can be calculated in the 
factorization approximation where the vector meson is coupled to the vacuum by one of the 
weak currents. We use the numerical value of a-r = —0.953 x 10^'' from Ref. |^ . 



B. The parity-conserving amplitudes 

A description of the physical nonleptonic decay amplitudes B ^ B' + n can also be 
performed by using a lowest-order chiral analysis. Employing a chiral lagrangian truncated 
at lowest order in the energy expansion for the PV (or s-wave) amplitudes yields results 
identical to those discussed above for pseudoscalar mesons. However if one defines the lowest- 
order chiral lagrangian for PC (or p-wave) amplitudes, one finds that such an operator has 
to vanish since it has the wrong transformation properties under CP. Thus, the only allowed 
chiral lagrangian at lowest order can generate PV but not PC terms. 

The standard method to compute the PC amplitudes is the so-called pole model. As 
shown in Ref. [^, this approach can be motivated by considering the transition amplitude 



for the nonleptonic emission of a meson 

{B'MMH^\B) = J d'xe''^^9{x'){B'\[dA{x),HU0)]\B) . (65) 
Inserting a complete set of intermediate states, {|^)}, one can show that 

{B'M,{q)\H^\B) = -J d^xe''^^B'\[A'^{x,0),HUO)]\B) - q^Mt , (66) 
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where 

{B'\Af{0)\n){n\HUO)\B) 



SiPn - pb' - q) 



PB-Pn 



+ c>{pb - Pn - q)— 



P%-<f - Pn 



(67) 



While the first term in Eq. ( pBD becomes the commutator introduced in Eq. (pUD, the 
second term contains contributions from the ground state baryons which are singular in 
the SU(3) soft meson limit. These pole terms become the leading contribution to the PC 
amplitudes. We note in passing that in principle, such baryon pole terms can also contribute 
to the PV amplitudes, however, more detailed studies showed that their magnitude is 
only several per cent of the leading current algebra contribution. 

We begin by computing the p-wave amplitude of the A Nti decay since here we can 
compare with experiment. The contribution to the PC weak vertex coming from the baryon 
pole diagrams shown in Figs, ^(a) and (b) are given by 

-Btt = S'nntt ^Nk + A'aett ^NT. (68) 

where Ana and Any. are weak baryon baryon transition amplitudes. These quantities 
can be determined via current algebra/PCAC as before 

\\m{n%\H,^\A) = -!(n|[F,^o, //pv]|A) = ^{n\H,,\A) 
lim(vrV|ifp.|S+) = Z!(p|[F5,^^^^]|s+) = . (69) 

Then assuming no momentum dependence for the baryon S-wave decay amplitude and ab- 
sorbing the i factor in the definitions of Ana and Any, we obtain 

Ana = i{n\Hp^\A) = 2F^{7i%\Hp^\A) = ~V2F^{7i-p\Hp^\A) = -4.32 x 10"^ MeV 

^iVE = ^(p|i/pc|S+) = V2F,(A|^pv|S+) = -4.35 X IQ-^ MeV. (70) 
v2 

For the rj contribution the PC AN?7 term, shown in Figs, ^(c) and (d), can be written 

as: 
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B„ — gt^f^ri Ana + 5'AAr) ^a^a (71) 

mt^ — mN niiq — mt^ 

while for the kaon (Figs, ^e) and (f)), the expressions are: 

m„ — rrj-A rrin — rrv^ 

= dANK ^AfA — A'nek Ant. (72) 

— niA rriN — 

mp — m J 

= 25(nsk ^ ^jvs , (73) 

niN - rriT 

where we have used the relations A^to = —Ant and Apj^+ = \/2Ant- 
The expressions for the vector mesons are similar: 

«p = 3nnp—^——Ana + gL, ^ Ant (74) 
mA ~ ttin m-N ~ ^t 

= alnuo ^A^A + glf,ui -4ArA (75) 

m\ — rriN ttln — m\ 

PCV V 1 1 

Ck* = 5'Ink* Ank - fi-NEK* Ant (76) 

— mA — ms 



D-J'- = 2gl,^, Ant. (77) 

"^Af - rriT 

and the tensor coupling constants /?p, C^f'^ and -D^?'^ are obtained from the previous 
expressions by replacing the strong vector couplings with the tensor ones. The numerical 
values of all these coupling constants can be found in Table 

We point out that some studies have included meson pole diagrams of the form shown 
in Fig. ^ The contribution of these diagrams would be given by 

5'aNK 5 T^KTT (78) 

mi — mi 

K 7r 

where the meson meson weak transition amplitude Aktt can also be related via PCAC 
to the experimental amplitude for K ^ -ktx decay, yielding A^^^ = —2.5 x 10~^ MeV^ [44 
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There is considerable uncertainty regarding the phase between the meson and the baryon 
poles which lead some studies to adjust it to better reproduce the data. It has been argued 
that the presence of these meson pole diagrams is important to fulfil the requirements 
of the so-called Feinberg-Kabir-Weinberg theorem in the nonleptonic decays. On the other 
hand, counting powers of energy in a chiral analysis, one finds that while the baryon pole 
terms are of order the meson poles enter at next order, along with higher-order chiral 
lagrangians. In general, we found these contributions to be very small and have therefore 
neglected them in the following. In principle, SU(6)^ can be used as well to relate the weak 
meson — >■ meson pseudoscalar transition amplitudes with those of the vector mesons. The 
results, given here for completeness, are: 

A --^A 

Ak'p — ^Ktt 



Ak*u. = -^^K. . (79) 



VI. RESULTS 
A. vr-Exchange 

We begin our discussion by presenting the results using only the OPE part of the weak 
AA^ — > NN interaction. On one hand we would expect this meson to adequately describe 
at the least the long-range part of the transition potential, while on the other hand its 
contribution has minimal uncertainties since the weak ANvr vertex is experimentally known. 
It is therefore a good starting point to assess the significance of form factors as well as initial 
and final state correlations before including the other mesons in the potential. The results of 
our calculations with OPE only are shown in Table |IV| where the nonmesonic decay rate of 

is given in units of the free lambda decay rate (Pa). The uncorrelated results (FREE) are 
compared with computations that include initial AN short-range correlations (SRC), form 
factors (FF), and final-state interactions (FSI) separately for the central (C), tensor (T) 
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(adding to a total parity- conserving (PC) contribution) and parity- violating (PV) potentials. 
The free central term is reduced dramatically by the initial SRC, however, most of the 
uncorrelated central potential contribution is in fact due to the 5-function in the transition 
potential which is completely eliminated by SRC. Without the (5-function, the central part 
is reduced by about a factor of two. Including SRC, FF, and FSI gives a negligible central 
decay rate. In contrast, the contribution of the tensor interaction is reduced only 10% by 
SRC and by 20-35% once FF and FSI are included as well. Therefore, the contribution of 
the central term amounts to less than 0.5% of the total vr-exchange rate. This behavior has 



been found and discussed by other authors as well [p!5| , p^ , p9| . On the other hand, our PV 
potential yields about 40% of the vr-exchange rate, at variance with older nuclear matter 
results that reported either a 15% |1^ or a negligible [1^ PV contribution to the rate. 



In previous papers p7| , p8| we have demonstrated the sensitivity of the calculated decay 
rates to the form of the initial SRC. In particular, we found that older calculations using a 
phenomenological NN correlation function [p!5| , p0| to simulate AN SRC in the initial state, 
rather than SRC based on a realistic YN meson-exchange potential as is done here, tend 
to overpredict the amount of initial correlations. We find similar results for the final-state 
interactions. As shown in Table |rV|, the PV and PC rates are reduced by about 35% and 
50%, respectively, when FSI are included via a correlation function based on a realistic NN 
potential (last two columns), rather than the 20% and 40% reduction obtained with the 
phenomenological NN SRC function of Eq. (|21|). It is comforting to see that the variation 
between different realistic NN interactions, such as the soft-core Nijmegen and a modern 
version of the Reid potential, plays essentially no role. This behavior can also be understood 
from Figs. ^ and ^ where we compare the different correlated wave functions for several 
channels and a relative momentum of kj. = 1.97 fm~^. The figures demonstrate that the 
phenomenological correlated wave function overestimates the realistic NN wave functions 
at intermediate distances (0.5 — 1.5 fm) while it underestimates them at short distances. 
As a result, the phenomenological approach overestimates the decay rate by about 20%. 
Including realistic final state correlations leads to an interference between central and tensor 
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transitions. For phenomenological FSI the central and tensor contributions to the total rate 
can be seen in Table to add incoherently. Once a realistic NN potential is used this 
incoherence is replaced by destructive interference. While this is only a small effect for the 
pion due to the small size of the central potential term this interference is more significant 
for the vector mesons where the central transition amplitude is large compared to the tensor 
term. 



The second quantity of interest displayed in Table |3 which is sensitive to the isospin 
structure of the transition amplitude is the neutron- to proton- induced ratio r„/rp. Notice- 
able is the smallness of the ratio which is due to the Pauli Principle that suppresses the final 
T = l, L = 2, S=l state with its antisymmetrization factor (1 — )-^+'^+^. That excludes the 
tensor transition in the neutron- induced rate, which gives rise to nn (T = 1) pairs, and it is 
precisely this tensor piece which constitutes the largest part of the OPE diagram. However, 
this argument holds only for relative AN S-states. For nucleons in the p-shell there exits a 
relative AN P-state which contributes a small but nonzero amount of the tensor potential 
to the neutron-induced decay. Note that including initial SRC, FF and realistic FSI reduces 
the Tn/Tp further by about 40%. This is due to the elimination of the central potential 
for which we obtain a value for Tn/Tp of about 1/3. In principle, one would expect F^/Fp 
=1/2 for the central term due to the statistical factor of 1/2 that accounts for two identical 
particles in the final state. For this number becomes 1/2.4 since we have 5 neutrons 
and 6 protons. The remaining difference comes from different ^Sq (T = 1) and '^S'l (T = 0) 
final state wave functions which enter the various spin-isospin channels and, therefore, lead 
to slightly different An nn and Ap np transition amplitudes. 

The suppression of the central potential term due to SRC, FF and FSI also explains the 
difference between the uncorrelated and the fully correlated ratio of PV to PC amplitudes, 
PV/PC, shown in Table |V|as well. More relevant than this ratio, however, is the asymmetry 
parameter, a\, defined in Eq. (y). This quantity, which measures the interference between 
the PC and the PV part of the amplitude, can be accessed experimentally in contrast to the 
PV/PC ratio which is merely of theoretical interest. We find that this asymmetry parameter 
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is only mildly sensitive to initial SRC and FF but changes by more than a factor of two when 
realistic FSI are included. This observable thus clearly demonstrates that for its accurate 
prediction the use of a realistic NN potential to describe the interactions in the final state 
is imperative. Below we will use realistic FSI generated with the Nijmegen potential for all 
results that include final state correlations. 

Comparing our results obtained here with older nuclear matter computations 
||l^ , [T5|JT7| , p!8[| we point out that the fully correlated total rate in nuclear matter is predicted 
to be in the range of 1.85 - 2.3, thus overpredicting our shell model calculations by more 
than a factor of two. When a Local Density Approximation (LDA) is performed p|,p!4| the 
rate reduces to F^^ = 1-45 Fa for ]^C. Although this value is further reduced when the same 
A wave function as in the present work is used, the LDA result is still larger by about 20- 
40%, depending on the choice of the Landau-Migdal parameter which measures the initial 
state correlations. The remaining difference may not be entirely surprising since the LDA is 
expected to work well for heavier nuclei and begins to break down for s- and p-shell nuclei. 
We also find that our PV potential yields about 40% of the vr-exchange rate, at variance 
with nuclear matter results that reported negligible PV rates ||T^ . 



We conclude our discussion of the OPE-only by assessing the role of the relative AN 
P-state contributions. In shell model calculations such terms naturally arise for nucleons in 
p-shell and higher orbitals when one transforms from shell model coordinates to the relative 
AN 2-body system. For s-shell nucleons, where one has an S-wave in both the relative and 
center-of-mass (CMS) motions, the transition amplitude gives a maximum contribution at 
the back-to-back kinematics, ki = — k2, yielding a total CMS momentum K = 0. For p- 
shell nucleons one would expect the contribution from the initial relative AN S-state to be 
suppressed compared to that of the relative P-state since the CMS harmonic oscillator wave 
function is then a P-state and thus zero at K = 0. Surprisingly, we find that after integrating 
over all kinematics with ki ^ — k2, this relative L=0 term contributes about 90 % to the p- 



shell rate |TP|. Thus, once the whole phase space is included, most of the total decay rate of 



p-shell nucleons still comes from the relative AN S-state. Furthermore, neglecting the relative 

30 



P-state contribution leaves the ratio r^/Fp unaltered while the asymmetry parameter oa is 
reduced by 10%. 

B. TT- and p- Exchange 

In this section we begin examining the role of additional mesons. As discussed above, 
one is faced with the immediate difficulty that none of the weak couplings involving heavier 
mesons can be accessed experimentally. Thus, one is required to resort to models which 
in this case involve considerable uncertainty. Table presents our results for the p-meson 
exchange alone as well as for the vr- and p-exchanges combined. Since both the ANvr- and 
the ANp-couplings are obtained within the same model there is no sign ambiguity. As 
noted before, the central potential can now be divided into a spin- independent (C) and a 
spin-dependent (SS) piece, which are shown separately. In contrast to the pion case, the 
factor in front of the Yukawa function in the SS central part of Eq. ( ^Bf ) enhances this 
contribution which then becomes comparable in magnitude to the piece containing the delta 
function. The two terms interfere destructively and yield a SS central part that is about half 
the size of the tensor contribution. In Ref. we noted that SRC reduce both the central 
C- and the SS-part of the p-meson contribution without 5-function by almost a factor of 10, 
compared to a factor of 2 in the vr case, reflecting the much shorter range of the p-exchange 
diagram. Similarly, the tensor interaction of the p is reduced by a factor of 2.5, compared 
to a 10% reduction in the vr case, as soon as SRC are included. The additional inclusion of 
FF and FSI further reduces both the central and tensor rates by substantial amounts. As 
is evident from Table 0, the final result of the central p contribution exceeds the p tensor 
term by almost a factor of two. Due to the strong destructive interference induced by the 
FSI between the tensor and the central part the total PC rate turns out to be smaller than 
either term alone. 

In terms of the combined vr and p contribution we find destructive interference between 
the two mesons for the PC rate but constructive interference for the PV decay mode. While 
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the TT-only PC rate is reduced by 16% when the p is added, the vr-only PV rate is enhanced 
by about 17% even though the p-only PV rate is very small. These counterbalancing inter- 
ferences lead to a combined tt + p total decay rate that is very similar to that of the tt alone. 
The neutron to proton induced ratio r„/rp, on the other hand, is slightly decreased. This 
could be, in principle, surprising since the tensor rate of the p-exchange is not as dominant 



as it is in the vr case. As noted in Ref. [0], this is due to an interference pattern of the 



central C and SS amplitudes which is destructive for the n-induced and constructive for 
the p-induced mechanisms. This yields a central rate which is basically p-induced. The 
strongest change can be seen in the asymmetry parameter a\ which is reduced by more 
than a factor of 2. This reduction can be traced to the above mentioned interference pattern 
between the PV and PC rates which are measured by this observable. 

The p-meson was the first of the heavier meson which was included in several earlier 
calculations. The first attempt was due to McKellar and Gibson [|15| in a nuclear matter 
framework which evaluated the weak ANp couplings using SU(6) and, alternatively, a fac- 
torization model. In their approach, which neglected the PV couplings, the phase between 
the TT and p amplitudes was not determined and their final results varied dramatically with 
their different models for the weak coupling constants. In a more recent calculation []TB 



Nardulli obtains the PC couplings in a pole model approach similar to ours. Besides the 
ground state baryon pole he includes the |* baryon resonance pole terms as well as K* pole 



contributions that appear in Ref. |[18| , |42[| but have been omitted here. The weak baryon- 
baryon transition amplitudes for the resonance poles are taken from a pole model analysis 
of hyperon PC pion decays which uses an F/D ratio of —1 for the weak baryon transition 
amplitudes and adjusts the overall coupling to the experimentally measured p-wave vr decay 
rates. As pointed out in Ref. the most serious problem with this fit is that it employs a 
K ^ IT weak transition amplitude in its K pole graphs that is about an order of magnitude 
larger than the strength extracted from the weak kaon decay mode K ^ Tin [Q. The 
K* pole contributions are calculated using a simplified factorization approach in which a 
number of terms are neglected The PV couplings of Ref. 0] are computed in a pole 
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model approach that includes baryon resonance poles with negative parity, belonging to the 
(70, 1~) multiplet of SU(6). In order to obtain the weak baryon transition amplitudes the 
experimental hyperon s-wave vr decays are used as input. Therefore, his approach for the 
PV weak ANp vertex is considerably different from the analysis used in this study. For the 
sake of comparison. Table [V| also lists the p term calculated with Nardulli's weak coupling 
constants. The PC transition potentials turn out to be very similar in magnitude while the 
PV rate is larger by more than a factor of 10 for the p alone. This increase enhances the 
vr + p total decay rate by about 25%, while the F^/Fp and oa are reduced by roughly the 
same amount. One should point out that the close agreement in the PC terms is fortuitous 
since the baryon resonance pole terms which are not present in our approach contribute 
about 30% to the weak ANp tensor coupling of Ref. [|16[. From this comparison it becomes 



obvious that there is considerable uncertainty in the determination of the weak vector meson 
vertices. 

C. K- and K*-Exchange 

The results of our calculations for the exchanges of the strange mesons, K and K*, are 
shown in Table The kaon is the lightest meson after the pion with a strong coupling 
constant (^ank of comparable magnitude to g^^^, but of opposite sign. The total decay rate 
for the K-only exchange diagram amounts to about 15% of the tt term which is the largest 
contribution among all the heavier mesons. As discussed further below the kaon therefore 
significantly interferes with the vr-only rate. 

In contrast to the vr and the p which are isovector mesons and the rj and uo which are 
isoscalar, the K and K* can contribute to both the T = and T = 1 weak AA^ NN 
transition potential and therefore have two independent couplings, Ck and Dk [see Eq. 
(p9|)]. Due to their isospin structure it was pointed out several years ago that including 



the kaon exchange has the potential to strongly influence the F„/Fp ratio. Using a simple 
schematic model that ignores the spin structure it was shown that the ratio F„/Fp could be 
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estimated with the expression 

^0 + ^ 



(80) 



where Aq and Ai are the isoscalar and isovector couphng constants, respectively, which in 
our case are given by 

Ao = ^ + D^ 

A. = ^. (81) 



With the PV values of Table p| , we obtain F^/Fp = 4.6 using Eq. (^), confirming the result 



quoted in Ref. ^J. However, we obtain F„/Fp = 0.23 when the PC coupling constants of 
Table |T| are used. Our complete result, which considers the spin structure and includes 
both PC and PV amplitudes, turns out to be F„/Fp = 0.26 as shown in Table This much 
smaller result suggests that one cannot draw conclusions about the ability of the strange 
mesons to drastically increase the neutron to proton ratio. Moreover, we have calculated 
this ratio using PV and PC amplitudes only but retaining the spin dependence and obtain, 
respectively, the values 1.69 and 0.03, far away from the estimates made above using Eq. 

As can be seen in Table the F„/Fp ratio for the kaon only is larger than the corre- 
sponding value for the pion by a factor of 2.5, while the asymmetry parameter obtained is 
very small. As discussed in the previous section our framework for the weak baryon-baryon- 
meson coupling constants assumes the validity of SU(3) (and SU(6) in the case of the vector 
mesons). No attention has been paid to the effects of SU(3) symmetry breaking which is 
known to be of the order of 30%. These effects have been addressed in a recent work by 



Savage and Springer |^6[ in the framework of next-to-leading order chiral pertubation theory 
(xPT). They point out that understanding the weak NNK vertex could elucidate a problem 
regarding the nonleptonic vr decay of free hyperons. While the PV (s-wave) amplitudes of 
these AS'=1 decays are adequately reproduced at tree level, the corresponding PC (p-wave) 
amplitudes cannot be well described using coupling constants from the s- waves as input. A 
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one-loop calculation of the leading SU(3) corrections performed in xPT, found that 



these loop corrections can change the tree level prediction of the p-wave amplitudes by a 
disturbing 100%, thus raising questions about the validity of xPT in this sector. As an alter- 



native it was suggested that large cancellations may occur between tree-level p-wave vr 
decay amplitudes which would magnify the SU(3) breaking effects. The one-loop corrections 
to the weak NNK vertex found in ref. on the other hand, modify the tree-level p-wave 



amplitudes by only up to 30%. If an experimental signature for these SU(3) corrections 
could be found in the nonmesonic decay it would provide insight into the applicability of 
xPT to these reactions. Table ^ shows the results of our calculations performed with the 
Savage-Springer weak NNK couplings. As expected, the kaon rates are roughly a factor 
of two smaller since the improved constants are reduced by about 30%. The values of the 
r^/Fp ratio and the asymmetry, on the other hand, are barely affected because all pieces of 
the transition amplitude are reduced by about the same amount. 

The K* vector meson is the heaviest meson exchanged in our weak AA^ NN transition 
potential. Nevertheless, due to its large weak NNK* and strong ANK* tensor couplings it is 
more important than either the p or the u. The central and tensor potential contributions 
are comparable in size, however, due again to the interference generated by the realistic 
FSI that mixes S and D states, the total PC rate turns out to be of the same magnitude. 
The larger PV coupling constant yields a PV rate which is significantly greater than the 
corresponding rates for the p and the u exchange contribution. The total K*-only decay rate 
is seen to be about half of the K-only rate but twice as large as the p- and cu-only rates. Due 
to the relative magnitude of the central and PC potential terms we find a F^/Fp ratio for 
the K* diagram which is 0.5, about five times larger than that for the vr-exchange ratio. This 
is a dramatic illustration of the difference in isospin structure between the various mesons. 
The K* diagram also exhibits more than twice the asymmetry parameter of the K-meson. 

Coincident ally, the K + K* total rate is only slightly larger than the K-only total rate 
even though the interference between the K and K* is important, as can be seen from the 
separate channel contributions. The PV rate is more than twice the K-only result and, as 
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a consequence, the asymmetry parameter is enhanced dramatically. The F^/Fp ratio also 
turns to be quite large, although it remains to be seen how the interference between all the 
mesons actually affects the final result for the observables. This is discussed in subsection 
E. 



D. r]- and LJ-Exchange 



Table |VI1| presents our results for the isoscalar mesons alone. The //-meson exchange 



contribution is by far the smallest of the different mesons included in our potential. This 
may come as no surprise since it is known that including the rj in phase shift fits of NN 
potentials influences the parameters only marginally. The main reason for the small size 
of the 77-exchange lies in the magnitude of the strong NNrj coupling constant, which is not 
well determined but is known to be much smaller than both (7^^^ and Qank- We have used 
the value of the Nijmegen potential which is gj^j^^i = 6.4 even though the NN phase shifts 
are very insensitive to this coupling. Recent data on t] photoproduction on the nucleon find 



a much reduced value of around g^NNr; = 1-4 [0. Should these conclusions be confirmed 
then clearly the 77 can safely be neglected in both the strong NN sector as well as the weak 



AA^ NN transition potential discussed here. As shown in Table [VI I| , the behavior of the 
rj contribution follows that of the vr term: a negligible central term, the largest piece coming 
from the tensor potential, and a PV rate about half the size of the tensor term. Note that 
in contrast to the pion there is no charge exchange term for the proton-induced decay, thus 
the F„/Fp for the r] only is larger than that for the pion by almost a factor of four. 

The cj-meson exchange provides an interesting contrast to the p contribution since for 
the strong NNp vertex the vector coupling is relatively small and the tensor coupling is large, 
while the reverse is true for the strong NNtu couplings. The weak PC and PV couplings (see 
Table PTID , on the other hand, are comparable in size. This pattern of strong couplings is 



reflected in the distribution of decay strength seen in Table |V11|. The largest contribution 



comes from the spin independent central term which is about twice as large as the corre- 
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spending term for the p. The u tensor term, however, is smaller than the p tensor potential 
by about a factor of seven. The interference of the various terms yield a total PC rate for 
the uj meson that is very similar in magnitude to the p but generates a much larger Fn/Fp. 
Similar to the p the PV rate is negligible due to the small size of the weak PV ANo; coupling 
constant. No other models are available for weak vertices involving this meson. 

E. The full weak one-meson-exchange potential 

In this subsection, we explore the effect of including all the mesons discussed before on 
the weak decay observables. For this purpose we show in Figs. and ^ the contribution of 
the different mesons to the integrand of Eq. (|18D for relevant transition channels. Fig. |^ 
displays the tensor transition ^Di (T = 0) of the PC amplitudes since it yields the 

most important contribution for pseudoscalar mesons and gives rise to important interference 
effects when the contributions of the mesons are added in pairs of identical isospin. As is 
evident from the figure, the pion-exchange contribution dominates, not only in magnitude 
but also in range; a consequence of the pion being the lightest meson. As expected, the 
kaon provides the second-largest contribution with a range somewhat less than that of the 
pion, followed by the heavier mesons with an even shorter range. Note that the contribution 
of each isospin-like pair [{n,p), (K,K*), (^7,0;)] interferes destructively, thus the large tensor 
contribution of each pseudoscalar meson is partially cancelled by that of its vector meson 
partner, an effect that can also be explicitly seen in Tables 0, and |VII| , discussed above. 
We have not shown the integrands of the central transitions since they are very small for 
the pseudoscalar mesons. 

Significant interferences are also observed for the integrands of the PV transitions ^5*1 — >■ 
^Pi (T = 0) and ^Si '^Pi (T = 1), shown in Fig. ^ Again, we find the pion to be 
dominant among the mesons in the T = transition, while the contribution of the other 
mesons play a more important role in the T = 1 channel. 

Below we discuss results for the different observables characteristic of the weak decay. 
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We compare the results obtained with the Nijmegen strong couphng constants with those 
obtained using the Jiihch strong couphngs given in Table Although in principle the 
strong couplings also affect the PC weak; vertices through the pole model, our goal here is 
to assess, for one particular model of weak couplings, the effect of using strong coupling 
constants from two different YN potentials which fit the hyperon-nucleon scattering data 
equally well. 



The results in Table [VII1| again demonstrate the significance of the short-range correla- 
tions and form factors in the nonmesonic decay. Adding the heavier mesons without form 
factors and SRC (column FREE) leads to a total rate that fluctuates significantly, with the 
additional mesons giving an appreciable contribution to the vr-exchange rate. This behavior 
is considerably suppressed by short-range effects, as shown in the second column. The rate 
is especially sensitive to the inclusion of the strange mesons. While including the p-meson 
has almost no effect the addition of kaon exchange reduces the total rate by almost 50% 
when the Nijmegen strong couplings are used. The reduction is mostly compensated by the 
addition of the K*, yielding a rate 15% below the pion-only decay rate. The situation is 
similar for the t] and u, their combined effect on the rate is negligible. Thus, with Nijmegen 
couplings adding the heavier mesons gives a reduction of only 15%. The situation is shghtly 
different when the Jiilich strong coupling constants are employed; their omission of the rj 
and their larger K* and u couplings lead to a total rate 15% larger than the pion-only 
rate. This indicates that the results are sensitive to the model used for the strong vertices, 
although both results are consistent with the present experimental values. This sensitivity 
to the strong coupling constants is unfortunate since it will certainly complicate the task 
of extracting weak couplings from this reaction. Improved YN potentials which narrow the 



range of the strong coupling constants are required to reduce this uncertainty. Table |V111 



also shows the results obtained when the NNK weak coupling constants derived in next-to- 
leading order xPT are used. Due to the smaller value of the coupling constants the 
effect of the K meson is reduced and thus the total rate is increased by about 10%. 

The results for the ratio of the neutron- to proton-induced partial rates r„/rp are shown 
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in Table |Xl The neutron- to proton-induced ratio is, as expected, quite sensitive to the 



isospin structure of the exchanged mesons. It has been known for a long time that pion 



exchange alone produces only a small ratio |jT5[. While the role of the p is limited it is 



again the inclusion of the two strange mesons that dramatically modifies this partial ratio. 
Including the K-exchange which interferes destructively with the pion amplitude in the 
neutron-induced channel (see, for example, the T = 1 PV transition amplitudes of Fig. H) 
leads to a reduction of the ratio by more than a factor of three. The K*, on the other hand, 
adds contructively. Again, an indication of this behavior can be seen in Figs. |^ and |^. 
In the T = 1 PV channel, relevant for the n-induced rate, the K and K* amplitudes have 
the same sign, whereas in both T = channels the interference between the two strange 
mesons is destructive and, as a consequence, the p-induced rate is lowered with respect to 
the n-induced rate. Using the Nijmegen strong couplings constants leads to a final ratio 
that is 34% smaller than the pion-only ratio, while using the Jiilich couplings leave this 
ratio unchanged, due again mostly to their larger K* and uj couplings. Employing the weak 
NNK couplings calculated with xPT we obtain an increase of the Fn/Fp ratio by 17% with 
Nijmegen couplings while the ratio remains unchanged for the Jiilich model. 

Even though it is not an observable. Table M also presents the ratio of PV to PC rates to 



aid in the comparison with other theoretical calculations. Again, adding the strange mesons 
produces the largest effect, especially when using the Jiilich strong couplings. The final 
PV/PC ratio is larger by more than a factor of two compared to the pion-only ratio. The 
results quoted in Ref. [1^ are of the order of 1 and, therefore, closer to our results obtained 



with the Jiilich model. 

The intrinsic asymmetry parameter, ca, shown in Table |I^ is also found to be very 
sensitive to the different mesons included in the model. This is the only observable which is 
changed dramatically by the inclusion of the p, reducing the pion-only value by more than 
a factor of two. Adding the other mesons increases oa, leading to a result about 30% larger 
than for vr-exchange alone in the case of the Nijmegen couplings and 50% larger for the 
Jiilich model. The effect of using the weak NNK couplings from xPT is very small for this 
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observable. 



F. Comparison with Experiment 

Our final results for various hypernuclei are presented in Table ^ We find an overall 
agreement between our results for the nonmesonic rate and the experimental values, espe- 
cially when the xPT weak couplings for the K meson are used, which yield somewhat larger 
rates. 

It has been the hope for many years that the inclusion of additional mesons would 
dramatically increase the ratio of neutron- to proton-induced rates. Here we find the opposite 
to be true. The final ratio greatly underestimates the newer central experimental values, 
although the large experimental error bars do not permit any definite conclusions at this 
time. On the other hand, the proton-induced rate which has errors of the same magnitude as 
the total rate is overpredicted by our calculations by up to a factor of two. It is the neutron- 
induced rate which has been very difficult to measure accurately. It is somewhat surprising 
that while both individual rates appear in disagreement with the data their sum conspires to 
a total rate which reproduces the measurements. Other mechanisms that have been explored 
to remedy this puzzle include quark-model calculations which yield a large violation of the 
AI =1/2 rule 1^1] , ^ , and the consideration of the 3N emission channel [ANN NNN) 
as a result of the pion being absorbed on correlated 2N pairs . A recent reanalysis |^2[ , 



which includes FSI of the three nucleons on their way out of the nucleus via a Monte Carlo 
simulation, shows that the 2N-induced channel further increases the experimental error bars 
and leads to an experimental value compatible with the predictions of the OPE model. 
However, the same reference shows that a comparison of the calculated proton spectrum 
with the experimental one favors values of Tn/Tp =2-3. It is therefore imperative, before 
speculating further about the deficiencies of the present models in reproducing this ratio, 
to carry out more precise experiments such as the measurement of the number of protons 
emitted per A decay, suggested in Ref. 
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Regarding the asymmetry parameter, comparison with experiment can only be made 
at the level of the measured proton asymmetry. As discussed in Section II, this quantity 
is determined as a product of the asymmetry parameter Ap, characteristic of the weak de- 
cay, and the polarization of the hypernucleus, Py, which must be determined theoretically. 
The energy resolution of the experiment measuring the decay of polarized ^^C produced 
in a (vr"^, K'^) reaction [l^ was 5 — 7 MeV which did not allow distinguishing between 



the first three 1~ states. Before the weak decay occurs, the two excited states decay elec- 
tromagnetically to the ground state. Therefore, in order to determine the polarization at 
this stage, one requires: i) the polarization of the ground and excited states, together with 
the corresponding formation cross sections, and ii) an attenuation coefficient to account for 
the loss of polarization in the transition of the excited states to the ground state. In Ref. 
||53|| , hypernuclear production cross sections and polarizations have been estimated for the 
(vr^, K^) reaction in the distorted wave impulse approximation with configuration- mixed 
wave functions. We note that the sum of the cross sections for the two excited 1~ states 
amounts to 40% relative to the ground state peak, which is consistent with the (31 ±8)% 
obtained in a fit to the Brookhaven ^^C(7r+, K^))^C spectrum |E3. Using the values of Ref. 



for the polarization and cross sections of the 1~ states in together with the spin 
depolarization formalism of Ref. we obtain Py = —0.19. This value, together with 
Ap = 0.151 (Nijmegen) or 0.175 (Jiilich), determined from a\ using Eq. @, leads to an 
asymmetry A = —0.029 (Nijmegen) or —0.033 (Jiilich), which lies within the uncertainties 
of the experimental result. 

Using the same model for ^^B, which predicts equal formation cross sections for the 



2i and 22 states, we obtain a polarization of Py = —0.29. However, hypernuclear structure 
calculations by Auerbach et al |^6[ predicted strong configuration mixing which reduced the 
cross section of the lower 2"^ state by a factor of three relative to the higher one. This 



prediction was verified by a reanalysis of older emulsion data Taking these relative 

weights into account, we obtain the value Py = —0.43, which is the one used in Table |^ 
and leads to better agreement with the experimental asymmetry. Just as in the case of the 
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proton- to neutron-induced ratio, the present level of uncertainty in the experiment does 
not yet permit using the asymmetry as an observable that differentiates between different 
models for the weak decay. 

In order to avoid the need for theoretical input and access Ap directly, a new experiment 
at KEK [BB| is measuring the decay of polarized \iie, extracting both the pion asymmetry 
from the mesonic channel, A-^-, and the proton asymmetry from the nonmesonic decay, A. 
The asymmetry parameter a^- of the pionic channel has been estimated to be very similar 
to that of the free A decay [^, and, therefore, the hypernuclear polarization can now be 



obtained from the relation Py = Aj^-faT^-. This in turn can then be used as input, together 
with the measured value of A, to determine the asymmetry parameter for the nonmesonic 
decay from the equality Ap = A/ Py. This experiment will not only allow a clean extraction 
of the nonmesonic asymmetry parameter but will also check theoretical model predictions 
for the amount of hypernuclear polarization. 

Finally, we briefly compare our results to previous calculations. The only other shell 
model calculations we are aware of are those of Ref. and Ref. [^. The results of Ref. 



are in agreement with ours while the preliminary results reported by Dubach et al. 
appear to be very different. Their uncorrelated OPE-only rate for is listed as 3.4 which 
is about a factor of 2 larger than ours while adding initial SRC, FF, and FSI reduce this 
rate to 0.5. This amounts to a reduction factor of almost 7, in contrast to our suppression of 
roughly a factor of two. Furthermore, their correlated rate for ^He is listed as 0.9, almost a 
factor of two larger than the j^^C result. Unfortunately, no details are given in Ref. |jl8| that 
address these problems. We note, however, that there are some unexplained inconsistencies 
between their recent results of Ref. |jl8[ and what was reported ten years before in Ref. []50 



where the correlated vr-exchange decay rate for ^^C is 2.0 while the addition of the other 
mesons lowers this value to 1.2. In fact, these values are more consistent with their ^He 
results, as well as with the effect of short range correlations found in almost all studies of 
the nonmesonic weak decay either in nuclear matter or finite nuclei. Our results for F^/Fp 



again differ from what it is reported in Ref. |jT8|, where a value F„/Fp = 0.2 is obtained 
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for TT-exchange alone but 0.83 when all the mesons are included. However, their results 
in finite nuclei are, surprisingly, quite different from their nuclear matter results, namely 
Tn/rp = 0.06 for TT-exchange alone and 0.345 when all mesons are included. With regard to 
the asymmetry parameter, the nuclear matter results of Ref. [|I^ are qualitatively similar 
to our results for the Nijmegen couplings. They obtain a value = —0.192 for vr-exchange 
alone and —0.443 when all the mesons are considered. 

VII. CONCLUSIONS 

In this study we have presented calculations for the weak nonmesonic decay mode 
AA^ — > A^A^ of A-hypernuclei. In contrast to most previous investigations performed in 
nuclear matter this work analysed this hypernuclear decay in a nonrelativistic shell model ap- 
proach. The initial hypernuclear and final nuclear structure are taken into account through 
spectroscopic factors. All possible initial and final relative orbital angular momenta are 
included in the baryon-baryon system. Realistic initial and final short-range correlations 
(SRC) obtained from AN and NN interactions based on the Nijmegen baryon-baryon poten- 
tial are employed in order to treat the nuclear structure details with as few approximations 
and ambiguities as possible. Our calculations were performed in a one-boson-exchange 
model that includes not only the long-ranged pion but also contributions from the other 
pseudoscalar mesons, the rj and K, as well as the vector mesons p,uj and K*. The weak 
baryon-baryon- meson vertices were obtained using SU(6) and soft meson theorems for the 
PV vertices and the pole model for the PC vertices. The primary goal of this work was 
to reduce nuclear structure uncertainties as much as possible so that our framework can be 
used to extract these weak baryon-baryon-meson couplings. 

Total decay rates evaluated with the full weak OBE potential fall within 15% of the value 
obtained with pion exchange only and reproduce the experimental data. This is due to the 
interference between the contributions of the heavier mesons whose individual infiuence on 
the decay rate can be substantial. Including the kaon exchange alone reduces the total rate 
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by almost 50%, this reduction is compensated by adding the other mesons, specifically the 
K*. In contrast to previous studies we found little influence from the />-meson even when 
we used a different model for the weak ANp couplings, similarly the uj contributes at the 
10% level. The dominant contribution beyond the pion-exchange mechanism is clearly the 
kaon exchange, followed by the K* which tends to partially cancel the effects of the kaon. 
It is therefore imperative that future studies include both strange mesons simultaneously. 
The importance of kaon exchange makes it possible to see the effects of modifying the weak 
NNK couplings by loop contributions from next-to-leading order xPT. Including these loop 
graphs leads to a reduction of the NNK couplings from their tree-level value up to 50%, 
which in turns modifles the rates by up to 20%. Future experiments should be able to verify 
this effect. 

We found the dominance of strange mesons to be even more pronounced in the partial 
rates and their ratio. Including the kaon reduces this ratio by more than a factor of three, 
which again is compensated by the K*. Furthermore, this ratio turns out to be sensitive 
to the choice of strong coupling constants as well. Using the Nijmegen strong couplings 
reduces the ratio by 30% from its pion-only value while the use of the Jiilich strong couplings 
leads to a change of only a few percent. This flnding indicates the need for improved YN 
potentials with better determined strong couplings at the hyperon-nucleon-meson vertices. 
Both theoretical values are far away from the experimental data, even though the error bars 
are still large. It appears to be impossible to reconcile these discrepancies within a one- 
boson-exchange potential. In order to approach the experimental values the weak couplings 
of the heavier mesons would have to be unreasonably large which would yield very large 
total rates incompatible with the data. If future experiments with improved partial rates 
confirm the present trend new mechanisms of a different kind would have to be introduced to 
resolve this puzzle. In contrast to the previous observables we found the proton asymmetry 
to be very sensitive to the p-exchange while the infiuence of the kaon is more moderate. This 
polarization observable is therefore an important addition to the set of observables since its 
sensitivities are different from the total and partial rates. The results of our calculations are 
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within the very large current error bounds. 



A. Criticisms 

Our study clearly indicates that further theoretical effort must be invested to under- 
stand the dynamics of the nonmesonic weak hypernuclear decay. Within the one-meson 
exchange picture it would be desirable to use weak coupling constants developed with more 
sophisticated approaches. A beginning has been made by Savage and Springer in their 
evaluation of the weak NNK couplings in next-to-leading order xPT and the effect has been 
found to be important. However, an understanding of the weak ANvr and SNvr couplings 
within the framework of chiral lagrangians is still missing. Furthermore, due to the im- 
portance of the K*-meson it would be desirable to recalculate its weak NNK* couplings in 



improved models as well. Several recent studies |^ have gone beyond the conventional 
picture of meson exchange and have developed mechanisms based purely on quark degrees 
of freedom. One should keep in mind, however, that such models have not always been able 
to reproduce the experimentally measured free hyperon decays. 

Another avenue that is currently being pursued is the validity of the AI=l/2 rule in the 
AA^ NN process. While this empirical rule is well established for the free hyperon and 
kaon decays there is some indication that it could be violated for the AA^ — > A^A^ process. 
Within the framework of SU(3) and soft meson theorems the weak vertices of NNK and 
AN?7 are related to the observable ANvr decay, therefore, one would expect small A/=3/2 
contributions for these mesons. On the other hand, the vector meson vertices can receive 
substantial contributions from factorization terms which have been shown not to fulfil the 



AI=l/2 rule [^. The attractive feature of these additional terms is their strong influence 
on the ratio Tn/Tp. 

On the level of implanting the basic AA^ — > A^A^ amplitude into the nucleus uncertainties 
have been minimized our study by treating each ingredient as well as possible. Neverthe- 
less, within the framework of the impulse approximation and the shell model the short 
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range correlations, the spectroscopic factors and the single particle wavefunctions still come 
from separate sources. This dilemma can be avoided in rigorous few-body calculations with 
realistic wave functions. The nonmesonic decay of the hypertriton can be calculated us- 
ing correlated three-body hypernuclear wave functions for the initial hypertriton state and 
continuum Faddeev solutions for the three-nucleon scattering state Thus, all nuclear 
structure input is generated from the same underlying YN- and NN-potentials, eliminating 
the ambiguities of the shell model approach. It is therefore of utmost importance to pursue 
experimental measurements of the nonmesonic decay of the hypertriton. 



B. Outlook 

On the experimental side, it is critical to obtain new high accuracy data soon. Improved 
partial rates for the proton- and neutron- induced decay modes are especially important. Of 
help would be to not only measure rates but also exclusive spectra of the decay products. 
Such distributions would be significant to disentagle the effects of the ANN — > NNN process 
from the two-body process discussed here. Beyond improving the present data base for the 
weak decay of A-hypernuclei, there are two more avenues which would aid our understanding 
of the weak AS=1 hadronic interaction. 

First, with the advent of new, high precision proton accelerators such as COSY in Jiilich, 
it may become possible to perform a direct study of the time-reversed process pn — > Ap |^ . 



While the very low cross sections present in this direct investigation of the AS=1 baryon- 
baryon interaction will be difficult to measure, high efficiency detection schemes should 
allow determining a branching ratio of 10^^^. Thus, the strangeness changing hadronic 
weak interaction could be studied similarly to the weak parity-violating NN interaction. 
The asymmetry of this reaction has been measured at several kinematics which are sensitive 
to different parts of the meson exchange potential. Furthermore, measuring the pn —>■ Ap 
process directly would give access to a number of polarization observables since the A is 
self- analyzing. 
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Secondly, the hypernuclear weak decay studies should be extended to double-A hyper- 
nuclei. Very few events involving these exotic objects — whose very existence would place 
stringent constraints on the existence of the elusive H-dibaryon — have been reported. 
Studying the weak decay of these objects would open the door to a number of new exotic 
A-induced decays: AA AN and AA T,N. Both of these decays would involve hyper- 
ons in the final state and should be distinguishable from the ordinary AN — > NN mode. 
Especially the AA — > AN channel would be intruiging since the dominant pion exchange 
is forbidden, thus this reaction would have to occur mostly through kaon exchange. One 
would therefore gain access to the AAK vertex. 

Even with the demise of KAON, the promising efforts at KEK with an improved measure- 
ment of the ^He decay, the continuing program at BNL, and the advent of the hypernuclear 
physics program (FINUDA) at DA$NE represent excellent opportunities to obtain new valu- 
able information that will shed hght onto the still unresolved problems of the weak decay of 
hypernuclei. 
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VIII. APPENDIX 



A. Coefficients {{V S)JMj\da\{LrSQ)JMj) 



Spin-Spin transition 



{{L'S)JMj\d^\{LrS^)JMj) = {2S{S + 1) - 3) 5l^l' Ss^s 



Tensor transition 



{iL'S)JMj\d^\iLrSo)JMj) = SIl' SsoS Ssi 



where the coefficients are given in Table ^ 



PV transition 



Pseudoscalar mesons 



{{L'S)JMj\da\{LrSo)JMj) = (-l)^+^-^'v^v^25o + 1^21,. + lV2STT 

(10L,.0|L'0) 



^ ^ 5*0 I L' Lf 1 



S 1 i 



2 J V 



So S J 



Vector mesons 



{{L'S)JMj\da\iLrSo)JMj) = i{-iy-^'+^QV6J2So + lJ2Lr + 1^25 + 1 



(10L^0|L'0) 



L' 1 ^ 
5*0 i/ ^ 



^ 1 1^ 



lie 

2 2'-' 

lie 

V 2 2 y 
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FIGURES 
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FIG. 1. Non-strange (a) and strange (b) meson exchange contribution to the AA'^ — > NN weak 
transition potential. The weak vertex is indicated by the circle. 
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FIG. 2. K- meson weak vertices for ppK'^ (a), pnK'^ (b), and nnK^ (c). 
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FIG. 3. Baryon pole diagrams contributing to the PC weak vertices in the AN 
transition ampUtude. 
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FIG. 4. Meson pole diagram contributing to the AN — NN transition amplitude. 
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CD'*0JO0J'*CD <D'*C\10C\1-*<D 

odoodod d d d ci d d d 
III III 

FIG. 5. ^Si Di coupled channel NN wave functions for a relative momentum of kr = 1.97 
fm~^, obtained with the Nijmegen93 (solid line) and the Reid93 (dashed line) interactions. The 
dotted line represents the phenomenological correlated wave function discussed in the text. 
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FIG. 8. Same as Fig. 7 for the PV Pi (T = 0) and ^Si Pi (T = 1) transition 
amplitudes. 
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TABLES 



TABLE I. Possible '^^'^^Lj channels obtained in the weak decay of p-shell hypernuclei 





Weak decay channel 


L' 


Strong FSI 


L 




(AN) 


(NN) 


(NN) 




Central 


'So 




'So 




'So 






'Si 




'Si, ^Di 


'Pi 




'Pi 




'Pi 


'Po 




'Po 




'Po 


'Pi 




'Pi 




'Pi 


'P2 




'P2 




'P2, 'F2 


Tensor 






'Di 




'Di, ^Si 


'Po 




'Po 




'Po 


'Pi 




'Pi 




'Pi 


'P2 




'P2, 'F2 




'P2, 'F2 


P.V. 


'So 




'Po 




'Po 


's. 




'Pi 




'Pi 


'Si 




'Pi 




'Pi 


'Pi 




'Su 




'Si, ^Di 


'Po 




'So 




'So 


'Pi 




'Si, ^Di 




'Si, ^Di 


'P2 




'D2 




'D2 


'P2 




'D2 




'D2 
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TABLE II. Constants appearing in weak transition potential for the different mesons (in units 
of Gpm^^ ^ 2.21 X 10-^) 



K. 



ii) 
C 



K 



(i) 
SS 



K, 



K 



PV 



TT 



2M 



2M 



A 



^NNTT 
■ 2M 



2M 2M 



2M 2M 



A 



' 2M 



K 



1 5ank 
2M 2M 



1 ^ANK 

2M 2M 



Sank 
2M 



2M 



2M 



2M 



2M 



9nNP "I" ^NNp 

' 2M 



V 



2M 



2M 



^NNoj "I" S'nnu; 



2M 



2M 



9nN(jJ ~I~ ffNNCc' 

2M 



K* 



S'ank* 



) ffANK* ~^ ffANK* 

'2M 2M 



ffANK* ~^ ffANK* 

"2M 2M 



ffANK* ~^ ffANK* 

2M 
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TABLE III. Nijmegen (Jiilich) strong and weak coupling constants and cutoff parameters for 
the different mesons. Tlie weak couplings are in units of GprriT^'^ = 2.21 x 10^''. For the kaon and 



the /9- meson we also quote the weak couplings obtained by Ref. |46] and Ref. |16|, respectively. 



Meson 



Strong c.c. 



Weak c.c. 
PC PV 



5NN7r — 13.3 
5as. = 12.0 



B^=-7.15 



(7NNr, = 6.40(0.) 
5AAr? = -6.56(0.) 



Bri=-U.3 



A^=1.80 



K 



5ank 



-14.1(-13.5) 



CPc=-18.9 



<7nek = 4.28(3.55) 



=-14.0 m 

=3.20 11 



Cj7=0.76 

=0.40 |6| 
Dj7=2.09 

=1.50 |6| 



iNNp 



3.16(3.25) 



13.3(19.^ 



ap=— 3.50 



-3.39 161 



/3p=-6.11 



-7.11 161 



gl^p = 0(0) 

gl^^ = 11.2(16.0) 



ep=1.09 



=3.84 PI 



UJ 



9aau 



9aau 



■ 10.5(15.9) 
: 3.22(0) 

7.11(10.6) 

-4.04(-9.91) 



a, 



-3.69 
/3^=-8.04 



e^,= -1.33 



K* 



9ank* 



T 

9ank* 



5nek* 



-5.47(-5.63) 
-11.9(-18.4) 
-3.16(-3.25) 



a': 



PC,V_ 
K* 

PCT 



D 



PC,V_ 



-3.61 
-17.9 
-4.89 



C^y=-4.48 



L>py=o.6o 
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5^^, = 6.00(7.87) 



L'^?'^=9.30 



TABLE IV. TT exchange contribution to the AA^ NN decay rate of PC 



FREE SRC SRC+FF SRC+FF+FSI 

phenom. Nijm93 



Eq. (|1D 



Reid93 



c 


0.282 


3.4 X 10-3 


1.3 X 10-2 


4.2 X 10-3 


3.3 X 10-3 


4.0 X 10-3 


T 


0.858 


0.781 


0.637 


0.685 


0.566 


0.579 


PC 


1.140 


0.785 


0.650 


0.689 


0.531 


0.547 


PV 


0.542 


0.447 


0.389 


0.421 


0.353 


0.345 


F/Fa 


1.682 


1.232 


1.038 


1.110 


0.885 


0.892 


Fn/Fp 


0.182 


0.113 


0.120 


0.118 


0.104 


0.100 


PV/PC 


0.476 


0.570 


0.598 


0.610 


0.665 


0.631 


OA 


-0.594 


-0.420 


-0.506 


-0.484 


-0.238 


-0.242 



64 



TABLE V. TT and p exchange contribution to the AA^ — > NN decay rate of j^^C. The values in 



parentheses have been calculated using the weak coupling constants by Nardulli [ p!^ . 





TT 


P 


■K + p 


C (C) 





0.020 (0.019) 


0.020 (0.019) 


C (SS) 


0.003 


0.014 (0.016) 


0.006 (0.007) 


C (Total) 


0.003 


0.045 (0.047) 


0.028 (0.030) 


T 


0.566 


0.027 (0.032) 


0.373 (0.358) 


PC 


0.531 


0.021 (0.023) 


0.445 (0.428) 


PV 


0.353 


0.008 (0.096) 


0.414 (0.635) 


r/PA 


0.885 


0.029 (0.120) 


0.859 (1.063) 


Pn/Pp 


0.104 


0.076 (0.097) 


0.095 (0.063) 


OA 


-0.238 


0.036 (0.046) 


-0.100 (-0.008) 
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TABLE VI. K and K* exchange contribution to the AA^ NN decay rate of ]^C. The values 



in parentheses have been calc 
next-to-leading order in xPT 


ulated using 
4(:1. 


the the NNK weak coupling 


constants obtained from 




K 


K* 


K + K* 


C (C) — 
C (SS) 0.004 (0.002) 
C (Total) 0.004 (0.002) 
T 0.083 (0.038) 
PC 0.093 (0.044) 
PV 0.040 (0.018) 


0.019 
0.092 
0.038 
0.038 
0.037 
0.023 


0.019 (0.019) 
0.130 (0.122) 
0.063 (0.058) 
0.015 (0.005) 
0.082 (0.050) 
0.091 (0.061) 


P/Pa 0.133 (0.062) 


0.060 


0.173 (0.111) 


Tn/Tp 0.263 (0.272) 


0.500 


0.647 (0.760) 


OA -0.080 (-0.090) 


-0.192 


-0.426 (-0.532) 


TABLE VIL rj and uj 


exchange contribution to the AA^ — > A^A^ decay rate of a^C. 




n 




ri + uj 


C (C) 
C (SS) 
C (Total) 
T 

PC 
PV 


0.001 
0.001 
0.005 
0.006 
0.003 


0.045 
0.009 
0.036 
0.004 
0.024 
0.002 


0.045 
0.016 

0.036 
1.5 X 10""^ 
0.035 
0.005 


P/Pa 


0.009 


0.026 


0.041 


Pn/Pp 


0.383 


0.235 


0.183 


OA 


-0.114 


-0.086 


-0.134 
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TABLE VIII. Free and fully correlated nonmesonic decay rate of in units of the free A decay 
rate T\. The values in parentheses have been calculated using the Jiilich-B coupling constants at 
the strong vertex. 

FREE SRC+FF+FSI 

TT 1.682 (1.682) 0.885 (0.885) 

+p 2.055 (2.325) 0.859 (0.831) 

+K 1.336 (1.699) 0.497 (0.506) 

+K* 2.836 (3.821) 0.760 (0.902) 

+r] 2.467 (3.821) 0.683 (0.902) 

+UJ 2.301 (4.338) 0.753 (1.023) 



weak K-couplings 
from 



0.844 (1.104) 
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TABLE IX. Weak decay observables for ]^C. The values in parentheses have been calculated 
using the Jiilich-B coupling constants at the strong vertex. 



TT 0.104 (0.104) 

+p 0.095 (0.096) 

+K 0.030 (0.029) 

+K* 0.049 (0.070) 

+ri 0.058 (0.070) 

+UJ 0.068 (0.109) 



PV/PC a A 

0.665 (0.665) -0.238 (-0.238) 

0.930 (1.137) -0.100 (-0.052) 

2.413 (3.206) -0.138 (-0.074) 

1.797 (1.968) -0.182 (-0.202) 

2.249 (1.968) -0.200 (-0.202) 

2.077 (1.675) -0.316 (-0.368) 



weak 

NNK-couplings 0.080 (0.108) 1.678 (1.436) -0.302 (-0.350) 

from xPT [||] 
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TABLE X. 

1 1 1 j_ 1 

been caicuiated 
XPT [46|. 


Weak decay observables for 
using tne IN INK weak coupl 


various hypernuclei. The values in parentheses have 
ing constants obtained from next-to-leading order in 




iHe 


A B 


A ^ 


EXP: 


0.414 (0.467) 
0.41 ±0.14 


0.611 (0.686) 
0.95 ±0.13 ±0.04 H 


0.753 (0.844) 
1.14 ±0.2 
0.89 ±0.15 ±0.03 H 


Pn/Pp 

EXP: 


0.073 (0.089) 
0.93 ±0.55 


0.084 (0.099) 
1.04l°:i i 
2.16±0.58t[Jf i 
0.70 ± 0.3 [H 
0.52 ±0.16 m 


0.068 (0.080) 
1-331:^2^1 i 
1.87 ± 0.59l?;^2 § 
0.70 ±0.3 [51| 
0.52 ±0.16 g 


Pp/FA 
EXP: 


0.386 (0.428) 
0.21 ±0.07 


0.563 (0.624) 

o.3o;°:if i 


0.705 (0.782) 

o.3ilg:{? i 


OA 


-0.273 (-0.264) 


-0.391 (-0.378) 


-0.316 (-0.302) 


^(0°) 
EXP: 




-0.120 (-0.116) 
-0.20 ±0.10 


-0.030 (-0.029) 
-0.01 ±0.10 (10| 
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TABLE XI. Matrix elements of the tensor operator evaluated between generalized spherical 
harmonic states of definite J,L and S 



5 J ^, 



L' = J + 1 



Lr = J + 1 



-2(J+2) 
2J+1 



Lr = J-l 



6VJ(J+1) 
2J+1 



L' = J 



L' = J-1 



6VJ(J+1) 
2 J+1 



-2(J-1) 
2 J+1 
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